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We present results of lattice dynamics calculations of the P — V equation of state and the pressure 
dependence of the Raman-active E^g lattice phonon for p— H2 and o— D2 in a wide pressure range 
up to ~2 Mbar using our recently developed semi-empirical many-body potential, and density- 
functional theory. Comparison with existing body of experimental and theoretical results showed 
that the employed many-body potential is a reliable basis for high-precision calculations for phases 
I, II, and III of solid hydrogens. 
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An accurate determination of the equation of 
state (EOS) of solid hydrogens has been an im- 
portant research objective for decades. System- 
atic high-pressure studies were started in the sev- 
enties of the last century^-"- (see reviews^"- and 
references therein). At present these x-ray and 
neutron studies span the pressure range up to 
^2 Mbar-~— and temperature range up to 1000 
K. The highest compression reached in the EOS 
experiments is 10.4 for solid H 2 15 (7.6 for solid 
D^), essentially higher than that for solid he- 
lium (8.4)i£. The EOS data provide a funda- 
mental basis for examining intermolecular inter- 
actions, and for testing ab initio theories. A 
number of model intermolecular potentials'^ 
have been proposed based on the experimen- 
tal EOS data. Another experimental technique 
which complements x-ray and neutron diffrac- 
tion by providing direct information on inter- 
molecular interactions and vibrational dynamics 
is Raman scattering. The hep structure has a 
Raman-active optical mode (E 2g symmetry) in 
the phonon spectrum which corresponds to the 
out-of-phase shear motions in the two orthogo- 
nal directions in the ab plane. The frequency 
range of this Raman mode is extremely large, 
from 36 cm -1 at zero pressure^— to 1100 
cm -1 at 250 GPa. The Raman spectrum of 
solid molecular deuterium has been measured 
up to ^200 GPai2"— 2&£L. These measurements 
show that hep-based structures are stable in this 
pressure range. The calculations of the E 2g 
Raman frequency v{P) using various empirical 



potentials*^^ show that the result is highly 
sensitive to details of the potential used. There- 
fore, comparing the calculated and experimental 
v(P) is a hard test for any empirical potential (or 
for any other theoretical method, e.g. ab initio 
calculations) . It is essential that these properties 
are sensitive to different characteristics of the in- 
termolecular potential: while EOS is sensitive 
to the potential well depth, the Raman scatter- 
ing experiment probes the second derivative of 
the potential at the minimum. In our recent 
paper 30 we have proposed new semi-empirical 
isotropic potentials for H 2 and D 2 . Unlike the 
previous potentialsi2iLi&£L they include not only 
pair forces, but triple forces as well. The goal of 
the present paper is to perform detailed calcu- 
lations of the EOS and pressure dependence of 
Raman frequencies for H 2 and D 2 using our new 
potentials and to compare the results to available 
experimental data and theoretical results for a 
wide pressure range which spans the phases I, 
II, and III of solid p-H 2 and o-D 2 . 

As mentioned above, our potential includes 
pair (Up) and triple (U tr ) intermolecular forces 
(U tot = U p + Ut T ). This potential was de- 
signed in a manner similar to the potential for 
solid helium22^. It has the form of a sum of 
the pair Silver a- Goldman (SG) potential 17 (dis- 
carding the R~ 9 term) and three-body terms 
which include the long-range Axilrod- Teller dis- 
persive interaction and a short-range three-body 
exchange interaction. The latter was used in a 
Slater-Kirkwood form^^. Our potential also 



includes the translational-rotational interaction 
however we have found that its contributions 
both to EOS and Raman frequencies are neg- 
ligible. An explicit form and parameters of the 
employed potential are given in Ref. 30 . We re- 
strict ourselves to T = K, with the zero-point 
energy taken into account using the Einstein ap- 
proximation. A small pressure range (^ 0.5 
GPa) where quantum-crystal effects play a de- 
cisive role was excluded from consideration. 

The decomposition of the total energy into 
contributions from the pair forces (E p ), triple 
forces (E tr ), and the zero-point energy (E zp ) is 
presented in Fig. 1. The inset of Fig. 1 shows 
the respective decomposition of EOS. As can be 
seen, the interplay between these three contribu- 
tions is rather complicated and different for dif- 
ferent regions of the molar volume. At relatively 
small compressions Vq/V < 2, the ground-state 
energy (and consequently the EOS) is dominated 
by the zero-point energy E zp . The zero-point 
contribution decreases with rising compression, 
but it remains significant up to the highest 
pressures reached in EOS experiments. At the 
ten-fold compression the relative contribution 
E zp /E gs remains as high as 20%. It is interest- 
ing to note that for the eight-fold compression 
E zp /E gs coincides with that for solid 4 He, but 
for helium it increases with decreasing compres- 
sion more steeply and already at three-fold com- 
pression the zero-point contribution dominates 
in the ground-state energy^. The three-body 
attraction becomes appreciable at the compres- 
sions higher than two-fold which corresponds to 
pressures about 2 GPa. The relative contribu- 
tion of the tree-body forces \E tY \/E p monotoni- 
cally increases with pressure and for the eight- 
fold compression it reaches 0.5. 

There has been many attempts to propose 
effective pair potentials which would have the 
same softening effect as attractive many-body 
forces^^LL. To account for these effects Hem- 
ley et al modified the SG potentiate with a 
short-range correcting tern>^k This Hemley- 
Silver a- Goldman effective potential was shown 
to fit static compression data up to 40 GPa. 
The P(V) calculated with this effective pair po- 
tentials for higher pressures 14 increases far more 
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FIG. 1: Contributions of the pair E p and triple E tr forces, and 
zero-point energy E zp to the total ground-state energy E gs for 
solid p-H 2 . The inset shows the respective contributions to 
EOS. 



rapidly than in experiment. 

The calculated equations of state P(V) 
for solid hydrogen and deuterium are shown 
in Fig. 2 in comparison with DFT-GGA 
calculations^ and the experimental results from 
Refs.-^^^^^^k As can be seen, the semi- 
empirical calculation with the proposed many- 
body potential is in an excellent agreement with 
experiment in the pressure range 1 - 140 GPa 
(phases I and II). These results can be com- 
pared favorably with recently published EOS 
calculations^. From 140 GPa onwards, the the- 
oretical P(V) curve lies slightly below the ex- 
perimental one, and at the maximum pressure 
of 180 GPa (Phase III) the difference grows up 
to about 10%. The reason for this is the ne- 
glect of the higher order (n > 3) terms in the n- 
body expansion. The effect of the large-n terms 
increases with pressure, and at the metalliza- 
tion point the n-body expansion would converge 
extremely slowly. Methods based on the den- 
sity functional theory (DFT) within local den- 
sity approximation (LDA) and generalized gra- 
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FIG. 3: (Color online) Calculated and experimental Raman 
frequencies as a function of pressure for solid hydrogen and 
deuterium. Semi-empirical calculations for many-body po- 
tential: this work, p— H 2 (dashed line), calculations for the 
SG potential 5 ; o-D 2 (solid line). DFT-LDA theory (this 
work): H 2 o, D 2 •• Experiment (H 2 ): <M, x-^), o2S*2I ; (D 2 ): 



FIG. 2: (Color online) Calculated and experimental pressure- 
volume relations for solid H 2 and D 2 . Semi-empirical calcula- 
tions for many-body potential (this work): p— H 2 (solid line), 
o-D 2 (dashed line); DFT-GGA calculations (• • -)S&. Exper- 



iment (reduced to K): (H 2 : 



»); (D2: 



h— ; the inset shows the small-pressure range. Experiment (re- 



duced to K): (H 2 : 



(D 2 : C 3 , 



dient approximation (GGA) are somewhat of 
an opposite to the empirical potentials method. 
Indeed, the accuracy of the EOS from DFT- 
GGA— improves with the increase of pressure: 
in the pressure range 180-70 GPa the EOS from 
GGA practically coincides with the experimen- 
tal one and for P > 140 GPa the agreement is 
better than for our empirical potentials, but at 
lower pressures the ab initio results progressively 
underestimate the pressure, and GGA gives a 
strongly underestimated equilibrium volume of 
about 8 cm 3 /mol. The reason for this is twofold: 
first, GGA gives a poor description of the van 
der Waals forces, and second, DFT calculations 
ignore all quantum zero-point motions of nu- 
clei, including the distinction between ortho- and 
para-species. We also compare theoretical and 
the experimental values for the isotopic shift 
AP(V) = Pb 2 (V) - Pr 2 (V). In accordance 
with Ref.— we find that the empirical potentials 
strongly overestimate AP. 



The comparison between theoretical and ex- 
perimental pressure dependencies v{P) of the 
E 2g optical phonon Raman-active mode is pre- 
sented in Fig. 3. Since we could not find any 
DFT data of this mode in the literature, we 
have also calculated the E 2g Raman frequency 
of solid H 2 and D 2 in the Pca2i structure - one 
of the plausible candidates for the orientational 
structure of phases II and III - using DFT-LDA 
approximation. Our calculations were done us- 
ing the Full- Potential Linear Muffin-Tin Orbital 
(FP-LMTO) code RSPt^. 

Comparing the theoretical results with the ex- 
periment we see that similarly to that we had 
for the EOS at pressures lower than ^150 GPa 
the semi-empirical curves agree with experiment 
better than DFT calculations but at higher pres- 
sures the situation is reversed. The limiting pres- 
sures at which the semi-empirical approach still 
works are ^175 GPa while LDA has a fine agree- 
ment with the experiment for H 2 and with fre- 
quencies obtained for D 2 with the help of har- 
monic ratio of \[2 from 150 GPa up to the high- 
est considered pressures ~230 GPa. The fre- 
quencies calculated from the SG potential 17 de- 
viate from experiment even for very low pres- 
sures. The same is true 5 for the effective HSG 
pair potentiate!. Thus we have shown that 



while effective pair potentials work reasonably 
well for EOS up to 40 GPa, they fail for the dy- 
namical properties like Raman spectrum, where 
the explicit inclusion of the 3-body forces is nec- 
essary. 

In conclusion, we have calculated the EOS and 
the pressure dependence of the Raman-active 
E 2g mode using our recently proposed many- 
body potentials 30 , and compared the results to 
the experiment and previous semi-empirical and 
DFT calculations. Also, DFT-LDA calculations 
of the E 2g Raman frequency were performed. For 



phases I and II (P < 150 GPa) the proposed 
many-body potentials give excellent agreement 
with the experiment, much better then any pre- 
vious calculations. It proves that the new poten- 
tials are a reliable basis for high-precision calcu- 
lations of structure and dynamics of H 2 and D 2 
up to about 140 GPa. In particular, they provide 
a huge improvement over any effective two- body 
potentials, stressing the importance of includ- 
ing the 3-body forces. For the higher pressures 
(Phase III) the DFT approach is preferable. 
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